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Comment on PRB 59, 9195 (1999) and 62, 
14061 (2000) by D. S. Golubev and A. D. Zaikin. 

Golubev and Zaikin have presented calculationa^l which 
predict a finite dephasing time at zero temperature due 
to the electron-electron interaction. These calculations 
are in error because (i) they dajpot reproduce the re- 
sult of independent calculationsB'0, and {iillhe conclu- 
sion of ref. [ is-|»hysically inconsistentE^Q. In their 
recent pubhcationsElo Golubev and colaborators, (GZS) 
continue to insist that the Coulomb interaction between 
electrons in the regime of weak localization (WL) leads 
to a finite dephasing time T^piT) at zero temperature: 
Tip{Q) = tqzs- This contribution comes from the ultra- 
violet part of the electron-electron interaction and we 
will focus our discussion on the ultraviolet cutoff in the 
theory. 
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FIG. 1. The weak localization correction to the conduc- 
tivity is shown in Fig. 1(A). Figure 1(B) represents the ba- 
sic elements of the diagrammatic technique. The renormal- 
ized vertex is the solution to the diagrammatic equation, see 
Fig. 1(C). 

To our opinion, there is now overwhelming experme- 
ntal evidence against GZS statements. However, we be- 



lieve that since the debate is about a well-defined theo- 
retical problem, it has to be resolved by purely mathe- 
matical means without invoking the experimental results. 

We have already published a detailed discussion of the 
theoretical side of the problcmo. Taking into account 
the broad attention which the discussion attracted, and 
the current "antiperturbative" sentiments among some 
condensed matter theorists, we found it is worth writing 
this comment. We highlight one of numerous mistakes 
in their calculations to demonstrate that Refs. [ |l],||j|] 
presented an incorrect perturbative rather than a non- 
perturbative treatment of the problem. We address our 
comment to those experts in the field who still want to 
get a concrete answer to the question where GZS made 
the most crucial mistake. 

The weak localization correction to the conductivity 
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is determined by the diagram shown in Fig. 0, where the 
Cooperon C is defined in Fig. |[ Here Cd is the conduc- 
tivity of a d-dimensional sample and v is the density of 
states at the Fermi surface. Terms (a) - (f) in Fig. E de- 
scribe the effect of the electron-electron interaction on the 
Cooperon. The role of terms (a) in Fig. || is to cure the 
infrared divergence at small energy-momentum transfer 
{Q ^ l/^i/s) thus to restore the gauge invariancelj'Bfl. 
These terms do not diverge ultravioletly and thus are not 
important for our discussion. 

For remaining terms (b) - (f) we have 
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FIG. 2. The diagrammatic equation for the Cooperon C in the presence of electron-electron interaction. Note that the 
Cooperon and the vertex on r.h.s. of the equation contain the interaction as well. The basic diagram elements are defined 
in Figs. 0(B) and (C). The retarded, Keldysh and advanced fermionic propagator are labeled by R,A, and K respectively. In 
averaging procedure, one uses the relation G'^(e) = [G^(e) — G'*(£)] tanh(e/2r) for the exact Green functions of disordered 
system. 



The dephasing time is determined by r"-'^ = 
5](e, e,Q = 0), the Cooperon's self energy at £1 = £2 = £ 
and Q = 0. Expansion in Q leads to the correction to 
the diffusion constant, D, free of ultraviolet divergences. 



(See Ref. [ | for details). 



Taking only term (b) in Fig. into account. 



1 



S('^)(£,£,0) 



(3) 

we encounter an ultraviolet divergence and E('^)(e,e,0) 
reproduces tgzs from Eq. (71) of Ref. [ Q and Eq. (13) 
of Ref. [ ^. We used the following relation between 
the Keldysh component, L^{uj, q), and the retarded (ad- 
vanced) components L^(^^(ti;,q) of the screened inter- 
action propagator: L^(w,q) = coth(a;/2T)(L^((jj, q) — 
LA(^,q)). 

However, the contribution 



du! ImL^(cL',q) e — lo 

^ — T ~, 7~ "TT" LHiIlil . — . 

TT {2ttY (-iw + Dc(^) 2T 

(4) 



exactly cancels out this ultraviolet divergency. Thus, the 
self-energy E = Y.^^^ + E^'^-f) is determined by |w| < T 
(according to Eq. (0) e < TV 

Comparing Eqs. ^) and (|4|) with the result of Ref. 
we conclude the contribution of diagram (b) in Fig. 
included in the result of ref. [ |l|,||, (see Eq. (71) there), 
whereas the contribution of diagrams (c)-(f) in Fig. is 
omitted in all orders of perturbation theory in Ref. . 

In the framework of GZS's approach, the contributions 
from diagrams (a)-(b) can be associated with their 5/, 
and those from (c)-(f) with their iSr. The reason why 
GZ's results don't contain any contributions from (c)-(f) 
is because they simply neglect the contribution of iSj^ 
to the dephasing rate. They try to justify this (i) by 
claiming in Ref. [5,6] that Sr is purely real for all paths, 
and hence cannot contribute to dephasing, and (ii) by 
claiming in Ref. [1,5,6] that Sr^ — along classical time- 
reversed paths, and hence can be neglected in the expo- 
nent of their path integral when the latter is evaluated in 
the saddle-point approximation (in which only classical 
paths are considered). Finally, they claimed in Ref. [5,6] 
that (iii) they can reproduce perturbation theory by ex- 
panding the exponent of their path integral. 



All three claims are incorrect: (i) GZS obtain a real 
expression for Sr because in Ref. [1], because Eq.(43) of 
Ref.[l] neglects Poisson brackets that are needed when 
Fourier |-|transforming the pvV~ and V~ pv terms in 
Eq.(40)El. If the Poisson brackets are included, Sr con- 
tains an imaginary part (which in perturbation theory en- 
sures that the ultraviolet divergency in iSr cancels that 
of ^/). 

(ii) GZS obtain Sr ~ along classical time-reversed 
paths, because in Eq.(43) of Ref. they also [in addi- 
tion to neglecting Poisson brackets] , replacecl the density 
matrix by its " Wigner transform" , 
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This replacement is unjustified, as will be shown below, 
and leads them erroneously to conclude that Sr — 0. 

Claim (iii) is false, because after neglecting Poisson 
brackets, and using the "Wigner transform Eq. (||)", it 
is impossible to reproduce the perturbation result in any 
order, as it will be shown below. 

Let us now retrace the argument by which GZS claim 
in Ref. ^ that they can reproduce perturbation theory: 
GZS notice correctly that the diagrams (c)-(f) can be 
rewritten in coordinate representation as an amputated 
average [compare to Ref. [ ||, Eq. (A1,A2)] 



Im / drsdr4 
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X {Gf,{e)Gf,{e-uj)GUe)Gt2{e)). 



where the Cooperon legs corresponding to the average in 
Eq. (^) are cut, coefficient 27riyr^ is introduced to match 
with expressions (^) and (^), and G^4(e) is the Keldysh 
Green function 



G^4(e) = tanh 



2T 



(7) 



2 




GZS also present correctly [see Ref. I Eq. (A3)] 

Gf,ie) = J dvy {G^,, (e) - Gf,,) [6,,^ - 2p,,,] (8) 



where Si' 4 = S{rif — and 



(5i'4 - 2pi 



de 



Gf,,{e) 



(9) 



is related to the density matrix pi/4 of a non-interacting 
but disordered system. Substituting Eq. (^) into Eq. (||) 
and using the causality principle (analytical properties of 
the integral over cu) GZS obtain 

27rz/r2sir^^(£,£) =Ini J dr3dr4 J dry^Lf^iuj) (10) 

x(Gf3(e)G^i,(e - C.)[5P4 - 2piM]G?2(£)Gf2(e))amp. 

Equation (|l^) still contains all of the relevant contribu- 
tions, which, if they were calculated correctly (see below) 
would reproduce Eq. (^. It is important to emphasize, 
that the density matrix pin in Eq. ( p^ is a non-local op- 
erator constracted of Fermion operators taken at coinsid- 
ent moments of time but at different space points. For a 
disordered system, the density matrix p is a random long 
range object; the characteristic scale after which pi/4 can 
be considered to be local is jj^i — ^4! ~ \J D jT . Since the 
electron Green's functions G'^^^-' are also random and 
long range, the procedure of impurity averaging is not 
trivial, see below. 

We insist that GZS's uncontrollable approximation 
Eq.(||) of replacing the density matrix by its "Wigner 
transform", which they used to argue that Sr does not 
contribute to the dephasing rate, is impermissible. Once 
this approximation is used in Eq. (|l(]|), it is impossible 
to recover the correct result in any order of the pertur- 
bation theory: the "Wigner transform" neglects all the 



contributions to pi'4 from all electron trajectories except 
the straight line connecting points 1' and 4. On the other 
hand, those contributions are included in Eq. (|^). As the 
result of the frivolous replacement (||), GZS calculate two 
terms entering into one physical quantity (and canceling 
each other), namely iSr and Sj in their language, with 
different accuracy and arrive to their conclusions about 
zero temperature dephasing! It is neither a new physics, 
nor a non-perturbative treatment - it is just an incorrect 
evaluation of a pcrturbative contribution (10). 

The last point we want to clarify is how to obtain 
Eq. (I) from Eq. (|o|). A correct way is to use Eqs. (H) 
and (10) and perform the disorder averaging according to 
standard rules, see Fig. ^. A straightforward calculation 
(see Appendix) gives Eq. (|). Thus, the GZS's approach 
to introduce the density matrix into Eq. (^ instead of 
the Keldysh Green's function generates eleven diagrams 
out of original four. As a result of neglecting non-locality 
of the density matrix and substitution Eq. (||) all those 
diagrams are lost, i.e. scattering on all the impurities be- 
tween points 1' and 4 is not taken into account, compare 
Fig. I and Fig. 3 in Ref. [|. 
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APPENDIX 

In this appendix we consider diagrams, shown in Fig. 
I and demonstrate that Eq. leads to Eq. (g). We 
insert Eq. (^ into Eq. ( p^ to represent Eq. (^ in terms 
of the exact Green's functions, defined for a particular 
disorder realization. We have 

'\e,e)^lm / -^—^^ tanh — LI^(^, q) 



X $i2(q,W, £,£'), 



where 



(11) 
(12) 



$i2(q,w,e,e') = J dY^dr^drve''^'^''-'''^ x 

(Gf3(£)G^i,(e-a;)[G?4(£')-G^4(£')]G?2(£)G^2(£))amp. 

Following the standard technique for calculation of the 
disorder average $(q) we obtain diagrams, shown in the 
Fig. I 

The disorder averaged Green's function is 

G«'^(£, p + q) = , _^ ^ , , (13) 

£-^p -f VFq±z/2T 

where r is the mean free path time, — /2m — fi is the 
energy, corresponding to momentum p, /i is the chemi- 
cal potential and vp is the vector of the Fermi velocity 
directed along the momentum p. 



The diffuson T> and the Cooperon C for non-interacting 
system (shown by black boxes) are given by 



P(w,q) -C(c^,q) 
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and the renormalized vertex F (black triangular) is 

1 



(_Dq2 — iuj)T 



(14) 
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Now we are ready to write down the result for dia- 
grams in Fig. ^. For the Q = component of the Fourier 
transform 

<I>(q, Lo, e, e', Q) = y $i2(q, uj, e, e')e*^("^i-^^)d(ri - v^) 



we obtain 
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The sum of diagrams of /? and ^-type vanishes and for 
the sum of a and 7 type diagrams we obtain 



*(Q=o) 



— i{e — uj ~ e') 



Re 



1 



Dq'^ — i{e — e') 



— Attvit Re 



Dq^ - i{e - e'] 



(17) 
(18) 
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The first term in the last line is just a (5-function, which 
allows us to perform integration over e' . The second term 
vanishes after integration over lo due to analytical proper- 
ties of L^{uj) [causality principle]. As a result, we obtain 
Eq. i). 
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